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Abstract 



The lifetimes of the B~ and mesons are measured using the partially 
reconstructed semileptonic decays B —>■ D£~'PX, where D is either a or 
D*^ meson. The data were collected by the CDF detector at the Fcrmilab 
Tevatron collider during 1992-1995 and correspond to about 110 pb~^ oi pp 
collisions at s/s = 1.8 TeV. We measure decay lengths and extract the lifetimes 
to be t{B-) = 1.637 ± 0.058 ps and t{W) = 1.474 ± 0.039 1 {];[]^? ps, 
and the ratio of the lifetimes to be t{B~)/t(W) = 1.110 ± 0.056 where 
the first uncertainties are statistical and the second are systematic. 



PACS numbers: 13.20.He, 14.40.Nd 



I. INTRODUCTION 

Measurements of the lifetimes of the individual S-hadron species can probe their decay 
mechanism beyond the simple spectator model decay picture. In this model, all hadrons 
containing a heavy quark should have one identical lifetime, that of the quark. However, 
this picture does not hold in the case of charm hadrons; the lifetimes of and mesons 
differ by a factor of 2.5. Possible causes of lifetime differences include contributions from 
non-spectator decays, namely the annihilation and the VF-exchange processes, and so-called 
final-state Pauli interference effects. Obviously, these mechanisms play an important role 
in the decay of charm hadrons. However, they are expected to produce smaller lifetime 
differences between the B hadrons because of the larger mass of the b quark. 

In the past few years, the heavy quark expansion technique has been applied extensively to 
the calculations of inclusive decay rates of heavy hadrons, both spectator and non-spectator 
decays. It provides quantitative predictions for lifetime differences among the heavy hadrons. 
It is generally believed that there should exist a lifetime difference of order (5-10)% between 



6 



the B~ and 5° mesons. Bigi predicts that the B~ meson hfetime should be longer 
than the meson lifetime. However, Neubert and Sachrajda state that the sign of the 
deviation from unity cannot be predicted reliably. A much smaller difference, of order 1%, 
is predicted for the 5° and B^ meson lifetimes. 

Several direct measurements of B~ and B^ meson lifetimes have been performed recently 
by the e^e~ experiments and by CDF Indirect information has been obtained 

through the measurement of branching fractions [§. The precision of current measurements 
now approaches the level where the predicted small differences could be discerned, and im- 
provements in these measurements will provide a strong test of i?-hadron decay mechanisms. 

In this Article we report a measurement of the B~ and B^ meson lifetimes using partially 
reconstructed semileptonic decays. The data used in this analysis were collected in 1992-95 
with the CDF detector at the Fermilab Tevatron proton-antiproton collider at a center-of- 
mass energy ^/s = 1.8 TeV and correspond to an integrated luminosity of about 110 pb^^. 

In order to identify semileptonic decays of B mesons, events with a lepton (e~ or ^~ , 
denoted by i^) associated with a or D*^ meson are selected. (Throughout this Article a 
reference to a particular charge state also implies its charge conjugate.) The £^D^ candidates 
consist mostly of B^ decays, and the i~D*~^ candidates consist mostly of B^ decays. The 

mesons are reconstructed using the decay mode K^n^. The D*^ decays are 

reconstructed using the decay mode D*^ D^n^, followed by D° — K^tt^, K^ti^ti^ti^ or 
K^n^Ti^. About 6000 such decays are reconstructed in the data sample. The decay length 
distributions are measured and the lifetimes are extracted after correcting for the relative 
admixtures of B~ and B^ mesons in the samples. The results presented here supersede a 
previous CDF measurement [Q, since the part of the data sample used here is the same as 
that of Ref. |]. 

II. CDF DETECTOR AND TRIGGER 

The CDF detector is described in detail elsewhere [0. We describe here only the detector 
components most relevant to this analysis. Inside the 1.4 T solenoid the silicon vertex de- 
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tector (SVX) and the central tracking chamber (CTC) provide the tracking and momentum 
analysis of charged particles. The CTC is a cylindrical drift chamber containing 84 measure- 
ment layers. It covers the pseudorapidity interval \r]\ < 1.1, where rj = — ln[tan(0/2)] [Q. The 
SVX consists of four layers of silicon micro-strip detectors located at radii between 3.0 and 
7.9 cm from the interaction point and provides spatial measurements in the r-(p plane with a 
resolution of 13 /im. It gives a track impact parameter resolution of about (13+40/pt) /wni 0, 
where px is the transverse momentum of the track measured in GeV/c. The geometric ac- 
ceptance of the SVX is ~ 60%, as it extends to ± 25 cm from the nominal interaction point, 
whereas the position of the primary interaction vertices has an rms width of ~ 30 cm along 
the beam [z) direction. The transverse profile of the Tevatron beam is circular and has an 
rms spread of ~ 35 /xm for the data taking period in 1992-93 and ~ 25 in 1994-95. The 
Pt resolution of the CTC combined with the SVX is (j{pt)/pt = [(0.0066)^ + (0.0009 pt)T^^- 
Electromagnetic (CEM) and hadronic (CHA) calorimeters with projective tower geometry 
are located outside the solenoid and cover the pseudorapidity region \ri\ < 1.1, with a segmen- 
tation of Aip = 15° and Arj ~ 0.11. A layer of proportional chambers (CES) is embedded 
near shower maximum in the CEM and provides a more precise measurement of electro- 
magnetic shower profiles and an additional measurement of pulse height. Also, a layer of 
proportional chambers (CPR) is installed between the solenoid and the CEM, and samples 
electromagnetic showers at about one radiation length. Two muon subsystems in the central 
rapidity region are used for muon identification. The central muon chambers (CMU) are 
located just behind the CHA calorimeter, and the central upgrade muon chambers (CMP) 
follow an additional 60 cm of steel. 

Events containing semileptonic B decays are collected using inclusive lepton triggers. 
CDF uses a three-level trigger system, where at the first two levels decisions are made with 
dedicated hardware. The information available at this stage includes energy deposits in the 
CEM and CHA calorimeters, high px tracks found in CTC by a track processor, and track 
segments found in the muon subsystems. The Ej- threshold for the principal single electron 
trigger is 9 (8) GeV for the data taking period in 1992-93 (94-95), where Et = E sin 6, and E 
is the energy measured in the CEM. In addition, a track is required in the CTC with pt > 7.5 
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GeV/c that points at the calorimeter tower in ip. For the 1994-95 data taking period the 



CES was added to the trigger system |T0| . The electron trigger requires the presence of pulse 
height in the CES corresponding to an electromagnetic shower of 4 GeV or above. Also, the 
if position of the shower is available with a segmentation of Ay^ = 2°, and the CTC track 
is required to point at the shower. The single muon trigger requires a track in the CTC, 
corresponding to a particle with px > 7.5 GeV/c, and track segments in both the CMU and 
CMP systems that match the CTC track within 7.5° in (p. At the third level of the trigger, 
the event selection is based on a version of off-line reconstruction programs optimized for 
speed. The lepton selection criteria used in level 3 are similar to those described in the next 
Section. 



III. RECONSTRUCTION OF SEMILEPTONIC DECAYS OF B MESONS 

The analysis starts with identification of leptons, e~ or fi^. If an event contains a good 
lepton candidate, we look for the charm meson or D*^ produced in the vicinity of the 
lepton candidate, to be consistent with the semileptonic decay signature B — > i~vDX. A 
proper correlation between the lepton charge and the charm flavor, £~ with not £^ with 

is required. 



A. Lepton Identification 

The identification of electrons makes use of information from both calorimeters and track- 
ing chambers. To be specific we require the following: 



Longitudinal profile consistent with an electron shower, i.e. small leakage energy in 
the CHA. 



Lateral shower profiles measured in the CEM [|TT| and the CES |jT2[ consistent with 
electron test beam data. 
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• Association of a high px track with the calorimeter shower based on position matching 
and energy-to-momentum ratio. 

• Pulse heights in the CES and CPR consistent with an electron. 

Photon conversion electrons, as well as Dalitz decays of 7r° mesons, are removed by looking 
for oppositely charged tracks that have small opening angles with the electron candidate. 

Muons are identified based on the geometrical matching between the track segments in 
the muon chambers and an extrapolated CTC track. We compute the of matching, 
where the uncertainty is dominated by multiple Coulomb scattering in the detector material. 
We require < 9 in the r-(p view (CMU and CMP) and < 12 in the r-z view (CMU). 



B. Charm meson reconstruction 

To identify the £~D^ candidates, we search for K^rr^ decays near the leptons, 

removing events that are consistent with the D*~^ D^tt^ decay chain. The K^n^ 
decay is reconstructed as follows. We first select oppositely charged pairs of particles using 
CTC tracks, where the kaon mass is assigned to the particle with the same charge as the 
lepton (called the "right sign" combination), as is the case in semileptonic B decays. The 
kaon (pion) candidate is then required to have momentum above 1.5 (0.5) GeV/c, and to 
be within a cone of radius AR = 0.6 (0.7) around the lepton in 7]-ip space, where AR = 



[AriY + (A(y9)2. To ensure accurate decay length measurement, each candidate track is 
required to be reconstructed in the SVX, as well as the lepton track. To reduce combinatorial 
background, we require the decay vertex of the candidate to be positively displaced 
along its flight direction in the transverse plane with respect to the position of the primary 
vertex. The primary vertex is approximated by the beam position P,|l3. To remove events 



consistent with the decay chain D*^ D^n^, we combine additional positive tracks with 
the candidate and compute the mass difference (Am) between the D^tt^ and the 
assigning the pion mass to the tracks. The Am resolution is measured to be 0.74 MeV/c^. 
We remove the candidate if any track exists that gives a Am value between 0.142 and 
0.148 GeV/c^. The resulting K^tt^ invariant mass spectrum is shown in Fig. [l|(a). We 
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fit a polynomial background and a Gaussian distribution to the spectrum and find a mass 
resolution of 11.3 MeV/c^. Also shown by the shaded histogram is the mass spectrum for 
the "wrong sign" {K^n^ with i~) combinations, where no significant signal is observed. We 
define the signal region to be in the mass range from 1.84 to 1.88 GeV/c^. The total number 
of events in the signal region is 5198, and the background fraction is estimated from the fit 
to be 0.53 ±0.02. 

To identify i~D*~^ candidates, we search for D*^ —>■ D^n^ decays using two fully recon- 
structed decay modes, D° K~tt^ and — > K^7r+7r+7r^, and one partially recon- 
structed mode, —* K^ti^ti^. For the D'^ K^tx^ and K^tx^ti^ modes, we apply the same 
momentum and cone requirements to the kaon and pion candidates as in the B —>■ i~D^X 
reconstruction. For the K^ti^ti^ii^ mode, the kaon (pion) candidate is required to 

have momentum above 1.2 (0.5) GeV/c, and to be within a cone of radius 0.65 (1.0) around 
the lepton candidate. Also, we require the decay vertex of the D° candidate to be positively 
displaced with respect to the primary vertex in the D'^ K^h^h^ti^ and K^tc^h^ modes. 
For the fully reconstructed modes, the candidate has to be in the mass ranges 1.83 to 
1.90 GeV/c^ and 1.84 to 1.88 GeV/c^, respectively. For the partially reconstructed mode, we 
require the mass of a K^n^ pair to be between 1.5 and 1.7 GeV/c^; we do not reconstruct 
the 7r° and in the subsequent analysis treat the K~tt^ pair as if it were a D^. For each 
mode, we reconstruct the D*^ meson by combining an additional track, assumed to have 
the pion mass, with the candidate, and computing the mass difference. Am, between 
the D^Tr~^ and candidates. Figures |I](b)-(d) show the Am distributions. In Fig. [^(d) the 
peak is broadened because of the missing vr*^ meson. Also shown by the shaded histograms 
are the spectra from the "wrong sign" low-energy pion (D^ti^) combinations. We define the 
signal region as follows. The two fully reconstructed modes use the Am range 0.144 to 0.147 
GeV/c^, and the K^tc^tc^ mode uses the range Am < 0.155 GeV/c^. The numbers of events 
in the signal regions are 935, 1166, and 2858, respectively. 

We estimate the numbers of combinatorial background events by using the shapes of the 
Am spectra of the wrong sign (D^n^) combinations and normalizing them to the number of 
events in the Am sideband. The estimated background fractions are 0.09 ± 0.01, 0.18 ± 0.02 
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and 0.37 ± 0.02, respectively. They are summarized in Table |. 

It is possible that real or D*^ mesons are accompanied by a hadron h~ that was 
misidentified as a lepton, and such events can be included in the above samples. The hadrons 
can be either the decay products of the same B hadron that produced the charm meson or 
the primary particles produced in pp —>■ bbX and pp —>■ ccX events. We investigate this 
possibility by studying the wrong sign combinations, with or D*^, which cannot 
originate from B meson decays. We see no evidence for signal in these combinations. Based 
on this study we estimate the contribution of the D^*^h~ pairs to our signal to be (1.2 1 1;2)%, 
where possible charge correlations between the charm meson and the hadrons are considered. 
We ignore this background, and treat it as a systematic uncertainty. 



IV. DECAY LENGTH MEASUREMENT AND MOMENTUM ESTIMATE 

A schematic representation of the B meson semileptonic decay topology is illustrated in 
Fig. 1^. The B meson decay vertex Vb is obtained by intersecting the trajectory of the lepton 
track with the flight path of the candidate. The B decay length Lb is defined as the 
displacement of Vb from the primary vertex Vp, measured in the plane perpendicular to the 
beam axis, and projected onto the transverse momentum vector of the lepton-D° system: 

_ {VB-Vp)-pf'''' 

• 

To measure a proper decay length of a i? meson decay, we need to know the momentum of 
the B meson. In semileptonic decays, the B meson momentum cannot be measured precisely 
because of the missing neutrino. We use p^ to estimate the B momentum for each event, 
which results in a corrected decay length defined as 

T I £~ D° 

X = LBmB/PT ■ 

We call it the 'pseudo-proper decay length'. The residual correction between p^ and p^ 
is performed during lifetime fits we shall describe later. 

A typical resolution on this decay length x due to vertex determination is 50 fim, in- 
cluding the contribution from the finite size of the primary vertex. For subsequent lifetime 
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measurements, we use only those events in which the resolutions on reconstructed decay 
lengths X are smaller than 0.05 cm. Also we require the proper decay length of the 
meson, measured from the B meson decay vertex to the decay vertex, to be in the range 
from —0.1 cm to 0.1 cm with its uncertainty smaller than 0.05 cm. These cuts reject poorly 
measured decays and reduce random track combinations. In addition, we limit ourselves to 
events with reconstructed decay lengths x in the range between —0.15 cm and 0.3 cm. These 
cuts have been applied already for the charm signals shown in Fig. |I|. 

As mentioned above, we have used the momentum of the lepton-Z^*^ system, to 
calculate the pseudo-proper decay length. However, we still need to account for the missing 
momentum to measure B meson lifetimes. We define the ratio K of the observed momentum 
to the true momentum as 

The K distribution is obtained from a Monte Carlo calculation. The ISAJET event gen- 



erator WM is used for the production of the h quark, where the shape of the pt spectrum 



is modified slightly to match the QCD calculation in the next-to-leading order |T^. The 



fragmentation model by Peterson and others ||16| is used. The CLEO event generator is 
used to describe B meson decays. In particular, the semileptonic decays adopt the model by 
Isgur and others (ISGW) [|l^. A typical K distribution thus obtained has an average value 
of 0.85 with an rms width of 0.11, and is approximately independent of the ^'^ in the 
range of interest, which is typically 15 to 25 GeV/c. It is also independent of the decay 
mode except for the partially reconstructed mode —>■ K^ti^tt^, which has a slightly lower 
mean value (about 0.80) because of the missing 7r° particle. Two K distributions are shown 
in Fig. ^. 

The lifetime is determined from a maximum likelihood fit to the observed pseudo-proper 
decay length distributions. The likelihood for the signal sample is given by 

i 

where Xi is the pseudo-proper decay length measured for event i, and the product is taken 
over observed events in the sample. The first term in the likelihood function represents a 
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B decay signal event, while the second term accounts for combinatorial background events 
whose fraction in the sample is /bg- The signal probability density function J^sici^) consists 
of an exponential decay function ^exp(— ^) defined for positive decay lengths, smeared 
with a normalized K distribution D{K) and a Gaussian distribution with width saf 

/K / Kt\ 
dK D{K) e{x) — exp G{x) , 
CT \ CT J 

where r is the B meson lifetime, c is the speed of light, 9{x) is the step function defined as 
0{x) = 1 for X > and 0{x) = for x < 0, and the symbol "(8)" denotes a convolution. G{x) 
is the Gaussian distribution given by 

1 ( x"^ \ 

G{x) = -1= exp 



and (7j is the estimated resolution on Xj. The scale factor s is introduced as a fit parameter 
and accounts for a possible incompleteness of our estimate of the decay length resolution. 
The integration over the momentum ratio K is approximated by a finite sum 



dKD{K) Y,D{Kj)l^K, 
j 

where the sum is taken over bin j of a histogrammed distribution D{Kj) with bin width 
AX. 

The pseudo-proper decay length distribution of combinatorial background events, 
J--bg{x)-i is measured using mass sideband events, assuming that they represent the combina- 
torial background events under signal mass peaks. The functional form of the distribution is 
parameterized empirically by a sum of a Gaussian distribution centered at zero, and positive 
and negative exponential tails smeared with a Gaussian distribution: 

+ (/+/A+) e{x) exp(-x/A+) (8) G{x) 
+ (/-/A_) e{-x) exp(+a;/A_) ® G{x). 

The shape of the background function (parameters f± and A±) and the resolution scale factor 
s, as well as the signal lifetime cr, are determined from a simultaneous fit to a signal sample 
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and a background sample. We use the combined likelihood L defined as £ = £sig -^bG) 
where £bg = V^k^^^G^^k) and the product is taken over event k in the background sample. 
The amount of combinatorial background /bg is also a parameter in the simultaneous fit. 
This parameter is constrained by adding a term = |(/bg ~ (/bg))^/<^bg negative 
log-likelihood —I = —InC. The average background fraction (/bg) and its uncertainty o"bg 
are estimated from the signal mass distributions (Table |). 

The background sample for the i~D^ mode is formed from the sidebands, defined by 
the mass ranges 1.74 to 1.79 and 1.94 to 1.99 GeV/c^. For the D*^ samples we use Am 
sidebands: we use the right sign (D°7r+) sideband 0.15 < Am < 0.19 GeV/c^ for the two 
fully reconstructed modes, and 0.16 < Am < 0.19 GeV/c^ for the K^-k^tc^ mode. 

The background samples are summarized in Table 

The pseudo-proper decay length distributions of the background samples are shown in 
Fig. ^, together with fit results. The background parameter values and the resolution scale s 
determined from the fit are listed in Table |IT1|. The corresponding decay length distributions 
of the signal samples are shown in Fig. |^. We find the lifetimes to be ct{B) = 489 ± 15, 462 ± 
18,472 ± 19 and 449 ± 14 fim for the four modes, where uncertainties are statistical only. 

As a check of the procedure, we measure the lifetime using the proper decay length 
measured from the secondary vertex Vb to the decay vertex. The proper decay length 
distributions are shown in Fig. ^, together with fit results. The lifetime numbers are sum- 



marized in Table |^ The result is in reasonably good agreement with the world average 
value of 124.4 ± 1.2 /xm Wm. 



V. B AND MESON LIFETIMES 

In order to extract the B~ and 5° meson lifetimes, we must take into account the 
fact that the i~D^ and i^D*^ samples are admixtures of the two B meson decays. The 
semileptonic decays can be expressed as B —>■ i~V'D, where D is a charm system whose 
charge is correlated with the B meson charge. If only the two lowest mass charm states, 
pseudoscalar (D) and vector (-D*), are produced, the i~^D*^ combination can arise only 
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from the meson decay. Similarly, the £~D^ combination comes only from B~ meson 
decays, provided that the from the D*~^ decay is excluded. However, it is known that the 
above two states do not saturate the total semileptonic decay rates. All data indicate that 
higher mass charm mesons, D** states, as well as non-resonant D^*^7r pairs, are responsible 
for the rest of the semileptonic decays. We do not distinguish resonant and non-resonant 
components, and refer to both of them as D**. 

These D** meson decays can dilute the charge correlation between the observed final 
states and the parent B meson. For example, the _D**o meson decays to D^*^^tc~ as well as 
[)(*)^Ti^ final states, resulting in misidentification of B~ meson decays as B^ D*^£~uX. 
Nevertheless, and i~D*^ combinations are dominated by B^ and B^ meson decays, 

respectively. As described below, the contamination of the wrong B meson species is only 
at 10-15% level. This enable us to extract the two B meson lifetimes. 



A. Sample composition 

We estimate the fraction of B~ decays g~ in the i~D^ and i^D*^ samples as follows. 
The production rates of charged and neutral B mesons and their semileptonic decay widths 
are assumed to be equal. We also assume the D** meson decays exclusively to D^*^7r via the 
strong interaction, thereby allowing us to determine the branching fractions, e.g. D^*^^tt'^ vs 
/}(*)07]-+^ using isospin symmetry. We consider four factors affecting the composition. First, 
the composition depends on the fraction (/**) of the D** mesons produced in semileptonic 
B decays, 

_ B(B j-uD**) _ ^ _ B(B £-uD) + B(B ^ j-uD*) 
~ B(B ^ i-uDX) ~ B(B i-uDX) 

The CLEO experiment measures the fraction of exclusive decays to the two lowest mass 



states to be 0.64 ± 0.10 ± 0.06 |g. Thus, we find that /** = 0.36 ±0.12. A few experiments 
have recently observed some D** modes [^, but the sum of exclusive modes still does not 
add up to the total semileptonic rate. Second, g~ depends on the relative abundance of 
various possible D** states, because some of them decay only to D*7t and others to Dtt, 
depending on the spin and parity. This relative abundance is not measured very well at 
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present. Changing the abundance is equivalent to changing the branching fractions for Z)*7r 
and Dn averaged over various D** states. We define a quantity Py as 

^ ~ B{D** D*n) + B{D** Dn) ' 
where B denotes a branching fraction. We assume the relative abundance predicted in 



Ref. [IBl], which corresponds to Py = 0.78. We also consider the extreme values Py = 0.0 
and 1.0. Third, the composition depends on the ratio of the B~ and S° meson lifetimes, 
because the number of i~D^*^ events is proportional to the semileptonic branching fraction, 
which is the product of the lifetime and the semileptonic partial width. Finally, the sample 
composition depends on the reconstruction efficiency of the low energy pion in the decay 
D*^ D^TT^. If we miss the pion and reconstruct the meson, the D*'^ decay is included 
in the i~D^ sample and the sample composition is altered. The efficiency is measured to be 
e(7r) = 0.93lo;2i by studying the rates of i~D*~^ events with respect to £~_D° events. 

We also have to take into account the differences in the reconstruction efficiencies for 
the B —>■ i~uD, D* and D** decay modes. We examine this effect by using the Monte 
Carlo events we have used to obtain the K distributions. The ISGW model was used for 
semileptonic decays. We find that the first two modes show very similar efficiencies, while 
the last mode has an efficiency that is lower by about a factor of two. 

The dependence of the B~ fraction g~ on the parameters /** and Py are illustrated 
in Figs. and |. We find that g~ = 0.85 for the £-D° sample and g- = O.lOl^o 
for the i^D*^ sample when the two lifetimes are identical. The central values correspond 
to the nominal choice of the parameters, /** = 0.36, Py = 0.78 and e(7r) = 0.93. The 
uncertainties reflect maximum changes in g~ when the three parameters are changed within 
their uncertainties, namely /** to 0.24 and 0.48, Py to 0.0 and 1.0, and e(7r) to 0.72 and 1.0. 

We also note that the momentum correction factors {K distributions) need to be modifled 
when the sample composition parameters are changed. The K distributions for the decay 
B — s> i~uD** have lower mean values because of additional missing particle(s), and changing 
the amount of D** decays results in changes in the K distributions. 

There are other physics processes that can produce the lepton-D*^*) signature. The largest 
background comes from the decay of the 5° meson, B^ i^uD**^, followed by D**^ — 
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The contribution of this process to the lepton-D^*^ signal is estimated to be about 
2%. Other processes such as B ^ t^v^D^*'>X followed by — > ^"z/^z/,-, and B D^D'^*^X 
followed by i~X, are suppressed severely because of branching fractions and kinematic 

requirements on leptons. We have ignored these backgrounds here. Therefore the fraction 
of B^ mesons is given hj = 1 — . We treat effects of the physics backgrounds as a 
systematic uncertainty. 



B. Lifetime fit 

We can now determine the B~ and -B° lifetimes with a combined fit of the and 
£-£)*+ samples. The likelihood is given by 

^= n |n[(l-/BG)-^SIG(x.) + /BG-^BG(x.)] n-^BG(a;fc)|, 
sample \ i k ) 

where the product is taken over event i in each signal sample, event k in each background 
sample, and over the ^~D^ and i^D*^ samples. For each signal sample, we use a two- 
component signal distribution function given by 

-^siG(a;) = ^"-^slG(a;) + (1 - 9')^iiGix), 

where ^sig(^) -^sici^) represent the B~ and B^ meson components, respectively. The 
dependence of g~ on the lifetime ratio is taken into account during lifetime fits. 

The resuh of the combined fit is ct{B') = 491 ± 17 /im, ct{W) = 442 ± 12 /xm, 
where the quoted uncertainties are statistical, and are correlated with each other with a 
coefficient of —0.308. From these numbers we calculate the ratio of the lifetimes to be 
t{B-)/t{W) = 1.110 ±0.056. 

The pseudo-proper decay length distributions of the i~D^ sample and the combined 
sample are shown in Figs. and 10. The results of the combined fit are superimposed. 



C. Systematic uncertainties 



The sample composition is a source of systematic uncertainty in the B meson lifetime 
determination. We change each one of the parameters /**, Py and e(7r) to another value 

18 



while keeping others at their nominal values, compute the sample composition g~ and fit 
for the two B meson lifetimes. The results are listed in Table ^ We interpret the observed 
changes as systematic uncertainties. 

Other sources of systematic uncertainties considered in this analysis are described below. 
They are summarized in Table 

We have estimated the decay length distributions from real data using mass sidebands, 
thus minimizing model dependence. However, the assumed functional form may not be fully 
adequate to describe the true shapes. Thus, we have considered an alternative parameteri- 
zation that includes additional exponential terms; this has turned out to give only minimal 
changes in the result. 

Physics and fake lepton background processes are studied by adding their simulated decay 
length distributions to the background function. 

Other sources of systematic uncertainties include our estimate of the decay length reso- 
lution and of the B meson momentum. We have introduced a resolution scale factor s and 
find a value of about 1.35. We change this factor to 1.0 or 1.7, fix it at the value and repeat 
the lifetime fitting procedure. We assign the observed changes as an uncertainty. The mo- 
mentum correction {K distribution) is subject to some uncertainty too, because it depends 
on the kinematics of B meson production and of semileptonic decays. An alternative pt 
spectral shape of the b quark production was considered, based on a comparison of lepton 
Pt shape in the real data and Monte Carlo events. A simple V — A decay model was tried 
in place of the ISGW model to describe semileptonic decays. We apply these changes and 
obtain new K distributions, and repeat the lifetime fits. The observed changes are listed as 
a systematic uncertainty. In addition, the K distributions are somewhat dependent on the 
lepton momentum and on the cuts used for electron identification. We assign uncertainties 
due to possible incompleteness in the treatment of these effects. As stated earlier, the mo- 
mentum correction depends on the assumed amount of B decays to D** mesons. This effect 
is already accounted for in the sample composition uncertainty. 

Also, we have applied a loose cut on the decay length in some modes, and it introduces 
a slight bias (about 2.5 fim) toward a longer lifetime. Here we quote the number without 
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correction to the final lifetimes and assign a systematic uncertainty. Finally, a possible 
residual misalignment of the SVX detector and the stability of the position of the Tevatron 
beam are considered. Some of these uncertainties are common to the two B mesons and 
cancel in the determination of the lifetime ratio. All these effects are combined in quadrature 
to give the total systematic uncertainty. 



VI. FINAL RESULTS AND CONCLUSION 

We have measured the lifetimes of the B~ and mesons using their partially recon- 
structed semileptonic decays B — > £~VD'^X and B i~VD*~^X. Our final results are 

r(5-) = 1.637 ±0.058 1[J:[}^^ ps, 
t{W) = 1A7A± 0.039 +^omi PS, 



t{B-)/t{B^) = 1.110 ±0.056: 



0.033 
0.030 > 



where the first uncertainties are statistical and the second are systematic. The result is 
consistent with other recent measurements We combine this measurement with the 

CDF measurement |^ using fully reconstructed decays, 

t{B~) = 1.68 ±0.07 ±0.02 ps, 
r(5°) = 1.58 ±0.09 ±0.02 ps, 
t{B-)/t{W) = 1.06 ±0.07 ±0.02, 

and derive the following CDF average: 

t{B-) = 1.661 ±0.052 ps, 
t{W) = 1.513 ±0.053 ps, 
t{B-)/t{W) = 1.091 ±0.050, 

where the uncertainties include both statistical and systematic effects. There exists a small 
(about 3 yum) correlation in systematic effects between the two measurements, such as due 
to detector alignment, and it is taken into account in combining the results. 
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The ratio of the two B meson hfetimes differs from unity by about 9%, or two standard 
deviations. This agrees with the small difference predicted by theory. The result is also 
consistent with the current world average value of 1.03 ± 0.05 The meson lifetime 
is consistent with the B^ meson lifetime p2| within the uncertainty. 
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TABLES 



B Mode 




mode 


mass range Am range Events Background fraction 








(GeV/c2) (GeV/c2) 




if- 


-7r+ 


1.84 - 1.88 Not D*+ 5198 0.526 ± 0.018 


e-D*+ 


K- 


-7r+ 


1.83 - 1.90 0.144 - 0.147 935 0.086 ± 0.011 


l-D*+ 


K- 




1.84 - 1.88 0.144 - 0.147 1166 0.183 ±0.015 


e-D*+ 


K- 


-tt+ttO 


1.50 - 1.70 < 0.155 2858 0.366 ± 0.016 


TABLE L 


Definition of signal samples, numbers of candidates and estimated background 










B Mode 




mode mass range Am range Events 








(GeV/c2) (GeV/c2) 






K-TT+ 


1.74 - 1.79, 1.94 - 1.99 Not D*+ 7200 


e-D*+ 






1.83 - 1.90 0.15 - 0.19 1769 








^TT- 1.84 - 1.88 0.15 - 0.19 5030 


i-D*+ 




K-7r+7r° 1.50 - 1.70 0.16 - 0.19 3809 




TABLE IL 


Definition of background samples and numbers of events. 


B Mode 




mode 


scale s /+ A+ (/xm) /_ A_ (/xm) 




K- 


7r+ 


1.38 ±0.03 0.404 ±0.008 531 ± 12 0.136 ±0.007 240 ±10 




K- 


7r+ 


1.32 ±0.07 0.487 ±0.017 498 ± 21 0.136 ±0.014 240 ±22 




K- 


7r"'"7r"'"7r~ 


1.38 ±0.03 0.328 ±0.011 362 ± 12 0.058 ± 0.008 216 ± 21 




K- 




1.39 ±0.04 0.536 ±0.011 612 ± 17 0.098 ± 0.008 274 ± 20 


TABLE IIL 


Background shapes obtained from a simultaneous fit of signal and background 


samples. 








B Mode 






mode ct{B) (//m) cr(I)°) (/xm) 








K-Tr+ 489 ± 15 128.0 ± 5.3 








K-Tr+ 462 ± 18 133.8 ± 5.6 


e-D*+ 






K-TT+TT+TT- 472 ± 19 125.3 ± 5.2 


e-D*+ 






K'TT+TT^ 449 ± 14 127.5 ± 5.0 



TABLE IV. B and meson lifetimes measured for individual decay modes. Quoted uncer- 
tainties are statistical only. 
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/** Pv e(7r) g- 


CT (/im) 


cor r el. 


t{B-) 


I u I U ^ 


B- 




rr 

coert. 




0.24 0.78 0.93 0.899 0.064 
36 78 93 851 105 
0.48 0.78 0.93 0.796 0.155 


491.3 zb 16.3 
491.0 ± 17.3 
492.0 ± 18.9 


448.0 ± 10.7 
442.2 ± 11.6 
434.2 ± 13.3 


-0.187 
—0.308 
-0.461 


1.097 ±0.049 
1 110 ± 056 
1.133 ±0.067 


0.36 0.00 0.93 0.806 0.000 
0.36 1.00 0.93 0.858 0.133 


491.0 ± 17.5 
491.3 ± 17.4 


448.2 ± 9.7 
440.2 ± 12.3 


-0.105 
-0.360 


1.096 ±0.048 
1.116 ±0.058 


0.36 0.78 0.72 0.790 0.105 
0.36 0.78 1.00 0.874 0.105 


494.5 ± 18.7 
489.8 ± 16.8 


441.7 ± 11.8 
442.3 ± 11.5 


-0.357 
-0.290 


1.120 ±0.060 
1.107 ±0.054 


TABLE V. B~ and lifetimes from a combined fit of l~D^ and t'D*^ samples under various 
sample composition conditions. Quoted uncertainties are statistical only and are correlated between 
B~ and B^ . Also listed are their calculated ratios. 


Source 


(/xm) 


Contribution to 
cr(50) 
(/xm) 


t{B~) 
t{BO) 


Sample composition 
fraction (/**) 

D** composition (Py) 

Low energy pion reconstruction 
Background treatment 
Decay length resolution 
Momentum estimate 

h quark -pT spectrum 

B decay model 

Momentum dependence 

Electron cuts 
Decay length cut 
Detector alignment 


+1 

-0 

+1 

-0 
+4 
-1 

±5 

+7 

-5 

±4 
±4 
±6 
±5 

+0 

-5 

±2 




+6 
-8 
+6 
-2 

±1 

±5 

+7 
-5 

±4 
±4 
±6 

±5 

+0 
-5 

±2 


+0.023 
-0.014 
+0.006 
-0.015 
+0.009 
-0.003 

±0.015 
±0.002 

±0.016 


Total 


±13 




+le) 

-15 


+0.033 
-0.030 


TABLE VI. A summary of systematic uncertainties in the B and i?'-* lifetime measurement. 
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FIG. 1. Charm signals reconstructed in the vicinity of leptons I . Four modes are shown: (a) 
K-TT+ {non-D*+), (b) D*+ D°7t+,D^ K-Tr+, (c) D*+ D°Tr+,D° K'tt+t^+tt- 
and (d) D*^ D^7r~^, — > K~Tr~^Tr^. Plot (a) shows the K~Tr~^ invariant mass spectra, and (b-d) 
show the Am distributions. Shaded histograms show wrong sign combinations, and in (a) they are 
scaled by 0.5 for display purposes. 
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FIG. 3. Distribution of the momentum ratio K (see text) for B ^ i vD^X, followed by 
— K~'K^ and K^tt+tt'^ decays obtained from a Monte Carlo calculation. 
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FIG. 4. Distributions of pseudo-proper decay lengths for lepton-Z? background samples (points). 
Also shown by the curve is the result of lifetime fits. Four decay modes arc shown: (a) 
B i-i?D°X, K-Tr+ (non-D*+), and B i-i?D*+X, D*+ D^Tr+, followed by (b) 

K-Tr+, (c) 1)0 ^ K-TT+TT+TT- and (d) L»° ^ R-tt+tt^. 
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Pseudo— Proper Decay Length (cm) 



FIG. 5. Distributions of pseudo-proper decay lengths for lepton-Z) signal samples (points). 
Also shown are the result of lifetime fits, signal (dashed curve) and background (dotted curve) 
contributions, and the sum of the two (solid curve). The four decay modes (a-d) are the same as 
in Fig. |. 
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FIG. 6. Distributions of the D proper decay lengths measured with respect to the B meson 
decay vertex (points). Also shown are the result of lifetime fits, signal (dashed curve) and back- 
ground (dotted curve) contributions, and the sum of the two (solid curve). The four decay modes 
(a-d) are the same as in Fig. ^. 
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FIG. 7. Fraction of mesons in lepton-D^*^ samples as a function of the D** meson 
fraction /** in semileptonic B decays. Vertical lines show the range of CLEO measurement [20|. 
The relative abundance of various D** mesons is fixed to Py = 0.78 (see text). Low energy pion 
reconstruction efficiency is fixed to 0.93 (solid curves), 0.72 and 1.0 (dotted curves). 
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FIG. 8. Fraction g oi B mesons in lepton-D'-*-' samples as a function of the average D** 
branching fraction B{D** — > D*7r) or Py. Vertical line corresponds to the prediction of the ISGW 
model |18|. The D** fraction (/**) is fixed to 0.36 (solid curves), 0.24 and 0.48 (dotted curves). 
Low energy pion reconstruction efficiency is fixed to 0.93. 
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FIG. 9. Pseudo-proper decay length distribution of the i~D^ candidates (points). CTirvcs 
show the result of the combined fit with l^D*'^ candidates: The B~ component (dashed curve), 
the component (dot-dashed curve), and the background component (dotted curve). 
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FIG. 10. Pseudo-proper decay length distribution of the i~D*^ candidates (points). The three 
decay modes are combined. Curves show the result of the combined fit with i^D^ candidates: 

The component (dot-dashed curve), the B~ component (dashed curve), and the background 

component (dotted curve). 
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